p-brane dynamics in (N + l)-dimensional FRW universes: a unified framework 
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We develop a velocity-dependent one-scale model describing p-brane dynamics in flat homoge- 
neous and isotropic backgrounds in a unified framework. We find the corresponding scaling laws in 
frictionless and friction dominated regimes considering both expanding and collapsing phases. 
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I. INTRODUCTION 

Cosmological inflation [l|, Q is a period of acceler- 
ated expansion in the early universe providing the most 
plausible solution to the flatness, horizon and magnetic 
monopole problems and explaining the origin of large- 
scale struture. In the string theory inspired brane infla- 
tionary scenario 0-0] > the inflaton may be identified with 
the distance between two branes or a brane and anti- 
brane whose interactions originate the inflaton potential. 
In the latter case (brane-anti-brane inflation) the infla- 
tionary epoch ends as a result of brane collision and anni- 
hilation and, in the process, lower-dimensional branes are 
produced by the Kibble mechanism and appear as topo- 
logical defects. Although the production of 1-branes (i.e. 
cosmic strings) is favoured, higher-dimensional p-branes 
may also be generated [8, 9]. As a consequence, inflation 
may offer a natural mechanism for the formation of p- 
brane defect networks evolving in a higher dimensional 
spacetime. 

The large-scale dynamics of cosmic string networks 
in (3 + l)-dimensional Friedmann- Robertson- Walker 
(FRW) backgrounds has been extensively studied using 
a Velocity-dependent One-Scale (VOS) Model [IM1, 
which accurately describes the cosmological evolution of 
the networks' root mean square (RMS) velocity and char- 
acteristic length. In [13j . this model was generalized 
to isotropic and anisotropic (N + l)-dimensional back- 
grounds. Furthermore, a similar model EMS devel- 
oped for domain wall networks in isotropic backgrounds 
was shown to successfully describe the results of high- 
resolution field theory numerical simulations [l?} • In this 
paper, we develop a more general model describing the 
dynamics of p-brane networks of arbitrary dimension in 
(N + l)-dimensional homogenous and isotropic universes 
in a single framework. 

The paper is organized as follows. In Sec. |TT]we obtain 
the p-brane equation of motion in a (N + l)-dimensional 
FRW background and apply it to the particular case of 
cosmic strings. In Sec. IIIII we derive the VOS equations 
of motion for p-brane networks. In Sec. IIVI we inves- 
tigate the different frictionless scaling regimes that arise 



in collapsing and expanding FRW universes and in Sec. 
IVl we study the friction dominated regimes. We then 
conclude in Sec. [VT1 



II. p-BRANE DYNAMICS 

Consider a p-brane whose thickness is much smaller 
than its curvature radii (zero-thickness limit), so that its 
world-history may be represented by 



(1) 



where Ua with v = 0,1,..., p are the coordinates 
parametrizing the (p + l)-dimensional world-volume 
swept by the p-brane, uq is a timelike parameter and 
are spacelike parameter. The p-brane equation of motion 
may be obtained by minimizing the action 



(2) 



with respect to variations of Here, = gapx°k% t> 
and a p is the p-brane mass per unit p-dimcnsional area. 

Prior knowledge of the trajectory described by the 
brane throughout its evolution allow us to define a real 
scalar field multiplet, 0, in the p-brane world-volume, 
described by the Lagrangian 



c=x-vm. 



(3) 



where X — ~<jf ^ d^cjf^ /2, a comma represents a partial 
derivative, and V(4> e ) is the potential (on the remainder 
of this section we shall omit the index e). The poten- 
tial, V, needs to have, at least, two degenerate minima 
in order to admit p-brane solutions which can be made 
arbitrarily thin by appropriate tuning. The initial condi- 
tions may be set up in such a way that the scalar field <f> 
defines a new p-brane whose velocity coincides with the 
velocity of the original p-brane. The new p-brane may be 
regarded as a domain wall (a (p— l)-brane in the (p+ 1)- 
dimensional spacetime spanned by the brane during it 
time evolution), whose dynamics is identical to that of 
the original p-brane. By varying the action, 
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S = / C^jd p+1 



(4) 
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with respect to 0, one finds the equation of motion 



where 



f 



fJ 



(5) 



Consider a flat FRW universe whose line element is 
given by 



ds 2 = a 2 (r/)(-dr] 2 + dx • dx) 



(6) 



where a(j]) is the scale factor, rj — dt/a is the conformal 
time, t is the physical time and x are comoving coordi- 
nates. 

One may identify the timelike coordinate with the con- 
formal time (uq = 77) , so that goo = a 2 (rj) . Moreover, the 
velocity of the brane may be taken to be orthogonal to 
the brane itself and, therefore, perpendicular to all spa- 
tial parameters of the brane (temporal-transverse gauge 
conditions) 



Q 



9ol = 9io 



Q 



(7) 



Consider a set of local spatial coordinates 
(til, •••) Up, itp+i) such that the brane is locally a 
coordinate surface for which Up+i is constant and it 
moves along this direction. The spatial coordinates 
(ui,...,u p ) may be chosen in such a way that they 
form an orthogonal set which locally parametrizes the 
brane and, for latter convenience, whose coordinate 
lines coincide with its principal directions of curvature. 
Since the p-brane is embedded in a flat FRW Universe, 
this coordinate system exists in the vicinity of any 
non-umbilic point of brane's surface [l8| . By introducing 
the scale factors of the coordinate system |x )U ,| = h^, 
the metric components in this coordinate system can be 
written as 



i 2 h? , if i = j , 



if 3 ¥= i j 



(8) 



so that g = —a 2l > p+2) h\...h 2 +l . Eq. flS)) then becomes 

4> + pH<j>-V 2 u <f> = -a 2 V,j > , (9) 
with % = a j a and where 




(10) 



is the laplacian for this set of coordinates. Here, we con- 
sidered the zero-thickness limit, neglecting the variation 
of the scalar field <f> on the directions tangent to the brane 
[lU (d(j)/dui = for i = 1, ...,p). We also taken h p+ i = 1 
so that dUp+\ is the infinitesimal arclength along the it p +i 
direction. 

Given Eq. © and using the method described in detail 
in (20|, one obtains the equation of motion 



(1-w 2 ) [(p + l)Hv- k] = 0, 



(11) 



(12) 



and is the comoving curvature vector along the prin- 
cipal direction of curvature m and v = uv is the brane 
velocity. Note that k is obtained by projecting the curva- 
ture vectors along the velocity direction. In doing so one 
is losing information about the acceleration of the brane 
perpendicular to v (normal acceleration) but the evolu- 
tion of v only depends on the tangential acceleration of 
the brane. 



A. Cosmic strings 

The world history of an infinitely thin cosmic string 
in a flat FRW universe can be represented by a two- 
dimensional world-sheet with x = x(?7, a), obeying the 
usual Goto-Nambu action. The equations of motion can 
be written as 



x + 2H(l-x 2 )x - e-^e" 1 ^)' 
e = -2He± 2 , 



with 



x • x' = . 



,'2 \ 2 



1 - X 



*2 



(13) 
(14) 



(15) 



where dots and primes are derivatives with respect to 77 
and tr, respectively. Consider 



x x 

v = -. y = 7-7 



(16) 



where v(rj,a) = |x|. The left hand side of Eq. (fT3| is 
given by 

x + 2H (l - x 2 ) x = vv + vir + 2U (l - v 2 ) vv , (17) 

where v is perpendicular to v. Moreover, the right hand 
side of Eq. p^|) gives 



x')' = 



1 d fy 



7 dl \7 



1, dv „ 
7^ al 



(18) 



where we have taken into account that 7 = (1 — v 2 ) -1 / 2 , 
dy/dl — k _L y and the fact that the physical length 
along the string is given by dl = \dx\ — |x'|cier. Hence- 
forth, the component of Eq. (Tj"3"|) parallel to v yields 



v + (1 - v 2 ) {2Hv - k) = 0, 



(19) 



where n — k • v, which is equivalent to Eq. (|lip in the 
particular case with p = 1 . 
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It is interesting to contrast this with the corresponding 
result for domain walls, in which case the coefficient of 
the cosmological damping term is N (rather than 2) and 
k = v-^^" 1 (rather than v-k). Hence, the differences 
between the macroscopic evolution of cosmic string and 
domain wall networks are expected to increase with TV. 



III. VOS MODEL 

Let us now consider the case of a network of p-branes in 
a (A + l)-dimcnsional FRW Universe. The RMS velocity, 
denoted by v = ■J (v 2 ), is defined as 



/ v 2 pdV 
JpdV 



(20) 



where p is the p-brane energy density and V is the phys- 
ical volume. An equivalent definition for 1-branes would 
be v 2 = J v 2 jdl/ J "fdl. The characteristic length, L, of 
the network is defined as 



(21) 



L N-p 



where p = V^ 1 J pdV is the average brane density. An al- 
ternative but less useful definition of the physical length- 
scale of the network would be 



(Ph) = 



L 



N-p 
ph 



(22) 



where (...) denotes a volume weighted average. This defi- 
nition has the advantage that L p h is sensitive only to the 
spatial profile of the network, and is therefore indepen- 
dent of the velocity distribution. However, since a direct 
relation between L p h, p and v does not exist, L is often 
a more useful variable than L p h- In any case, L and L p h 
are in general very similiar, except if the p-branes are 
ultra- relativistic. 

Multiplying Eq. (TTTT) by v, making the volume average 
and then dividing by v, one obtains 



where 



v + (1 - v 2 ) [(p + 1) Uv - R] = . 



(v (1 - v 2 ) k) _ Jv (1 - v 2 ) npdV 
v(l-v 2 ) ~ v(l-v 2 ) J pdV 



(23) 



(24) 



and we have made the assumption that (i> 4 ) = v 4 (see 
10]). Eq. ([2"B")) may also be written as 



dv 
~dl 



+ (1-- 2 ) 



(25) 



where i^ = (p + 1)H is the damping lengthscale, H = 
Ha is the Hubble parameter, k — RL/a is a dimension- 
less curvature parameter equivalent (for p = 1) to that 
of the original VOS model for cosmic strings [l(| EH- 



The frictional force originated by the interaction of the 
branes with ultrarelativistic particles may be included in 
Eq. (j25|) . by introducing an extra term in the damping 
lenghtscale, i^ 1 = (p + l)H + l~ Y . The friction length- 
scale, if, will be properly defined in in Sec. [V] 

Energy-momentum conservation in a FRW universe 
then implies that 



v) = o. 



(26) 



assuming that the p-brane network is statistically ho- 
mogenous and isotropic on large scales. Here, V = 
V^ 1 J VdV is the average brane pressure. The equation 
of state parameter of the brane gas is given by [21( : 



w 



V 
P 



N 



[{p+l)v 2 -p} 



(27) 



We may also include an extra term in Eqs. (|2"T|) and (|2l))) 
to account for the energy loss due to p-brane collapse and 
friction 



dp 
~dt 



Hp [D + (p + 1) v 2 ] 



cv 

T 



's 



(28) 



where c > is the energy- loss parameter and D = N — p. 
Consequently, the equation of motion for the character- 
istic lengthscale of the network yields 



dL ttt L _ 2 
— = HL H v 2 

dt Di d 



D 



(29) 



Eqs. (|2"5j) and (|2T))) constitute a unified VOS model for 
the dynamics of p-brane networks in {N + l)-dimensional 
FRW universes. In the following section, we shall obtain 
the corresponding scaling laws. The reader is referred to 
; 16] for a more detailed discussion of the more specific 
domain wall scenario with p = N — 1 . 



IV. FRICTIONLESS REGIMES 

Let us assume that the dynamics of the universe is 
driven by a fluid with w = constant =/= —1 so that a oc 
t£, where f3 = 2/(N(w + 1)) and U > is the time 
elapsed since the initial singularity (if dt* = dt) or the 
time remaining up to the final singularity (if dt* = —dt) 
at =0. We shall consider six different models labeled 
by M?, where s = ± depending on whether dt = ±dt* 
and i — 1,2 or 3, depending on whether p3 < 0, < 
f3 < 1 or (i > 1, respectively. The models M^~, and 
AI^ represent expanding solutions with t* — either 
at the the big-bang (Af^ and M3 ) or at the big rip (for 
Mi). The models , and M3" represent collapsing 
universes with t* = either at the the big-crunch (M^ 
and M3) or at an initial infinite density singularity with 
a* = 00 (for Mf). 
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A. Linear scaling solutions (if — oo) 

If the friction lengthscale becomes negligible compared 
to the Hubble radius then Eqs. (|25p and (f2T)f may have 
a linear attractor solution. This attractor solution corre- 
sponds to a linear scaling regime of the form 



and 



constant . 



(30) 



with 



k(k + c) 



p(l-0)D(p+l) 



I (l-p)kD 
p(k + c)(p+l) 



(31) 

The conditions < v < 1 and £ > are sufficient to 
show that models and M-f do not admit linear scal- 
ing solutions, which would require a negative energy loss 
parameter c. 

The RMS velocity, v, of maximally symmetric p-branes 
with a Sp-i ® W topology oscillating periodically in a 
Minkowski spacetime is given by 



< 



p - i + 1 



(32) 



with v^ ain — 1/2 and u^ax = p/ (P+ 1); which correspond, 
respectively, to i = p — 1 (only one of the principal cur- 
vatures is non-zero) and i = (for fully spherically sym- 
metrical p-branes). For /3 = 0, the curvature parameter 
k must be equal to zero for a linear scaling solution with 
v < 1 to be attained. The expansion (collapse) of the 
universe hinders (aids) the velocity of the branes, lead- 
ing to a smaller (larger) RMS velocity and a curvature 
parameter, k, larger (smaller) than zero. Therefore, one 
expects that 



< v < v max , 

Vmin < V < 1 



for Mi} 



(33) 



for Mf and M 3 ~ 



On the other hand the characteristic lengthscale of the 
network is necessarily constrained by causality and, as a 
consequence, L is expected to be bounded by the particle 
horizon at any given time. In the case of models , Mi± 
and Mg" this implies that 



L < d 



H 



dt' 



t «(*') 11-01 



(34) 



with ti — or ti = — oo (depending on whether s = 
— , respectively), or equivalently 

v 2 <{k + c)- 2 . 



(35) 



Some remarks regarding the energy-loss term are nec- 
essary. For simplicity, we shall consider the most triv- 
ial case of flat p-branes. A moving p-brane spans a q- 
dimensional surface (with q = p + 1) which has N — q 
degrees of freedom. Hence, if N < 2(N — q) then two 



flat q-dimensional surfaces do in general intersect but 
that is no longer true if N > 2(N — q) (or equivalently 
p < (N — l)/2). Hence, if the p-branes are thin (p-brane 
thickness much smaller than L), the energy-loss param- 
eter c may be much smaller than unity if p < (N — l)/2. 
Note, however, that the linear scaling solution in Eq. ([50)1 
may be attainable for j3 > 1 — p/N, even if c = with 



k' 2 



P{i-P)D(p + i) 



_ (l-f3)D 

v= ywm- (36) 



As a consequence, even if c <^ 1, the p-brane network 
may be able to attain a linear scaling regime in the Mi^ 
model. However, if c = 0, a linear scaling solution is no 
longer possible in a collapsing universe. 



B. Inflation and superinflation 

In the case of models M^ and M-f the expansion is fast 
enough (a > 0) to decelerate the branes and make the 
RMS velocity arbitrarily small. As a consequence, the 
inflationary models Al£ and M£~ undergo a stretching 



regime described by the scaling laws 

Lex a, v (x (Hay 1 (x a^ 1/p 



0. 



(37) 



If the p-branes are the dominant energy component of 
the universe then 



P = 



N(l + w b ) D + (p + l)v 2 



(38) 



In order to accelerate the universe one needs p > 1 (or 
equivalently Wf, < w c = (2 — N)/N) and, consequently 



v < 



D 



p+ 1 



(39) 



We may then conclude that only domain walls (D = N - 
p — 1) are able to drive an inflationary phase. 



C. Ultra-relativistic collapsing solution 

In the case of model Mi^ (which represents a collaps- 
ing universe with k < and < P < 1), the comoving 
Hubble radius, decreases with time. As a conse- 

quence, the curvature scale of the p-brane will necessarily 
become smaller than and will tend to be confor- 

mally contracted in physical coordinates whilst travelling 
at ultrarelativistic speeds with 



L p h oc a , 



7 oc a 



-(p+i) 



(40) 



where is defined by 7 = (1 — v 2 ) 1 / 2 and L ~ 7 1 / £) L 



-{N+l)/D 



ph oc 

Hence, at ultra-relativistic speeds L ceases 



to be an accurate measure of L 



P h- 
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V. FRICTION-DOMINATED REGIMES 

The interaction of p-branes with the ultrarelativistic 
particles in a radiation fluid results in a frictional damp- 
ing of the form 



dv 1 /I 2n 



(41) 



where If oc a p X p+1 N /p la d oc a p+2 is the friction length- 
scale, A oc a and p ra d oc a - ^ -1-1 ) are, respectively, the 
typical wavelength and the energy density of relativistic 
particles. 



A. Expanding universe 



at a higher rate than (for p + 1 > 1//3) as the uni- 
verse expands and eventually overcomes the character- 
istic length, L. This implies that the Kibble regime is 
necessarily transient. 



B. Collapsing universe 

In the case of a collapsing universe, the evolution of 
the network ends in a friction dominated era, with the 
p-branes coming to a standstill in comoving coordinates 
and then being conformally contracted with p oc a . 
The background temperature and density will eventu- 
ally approach those of the brane- forming phase transition 
and, as a consequence, the branes will dissolve into the 
high density background. 



If, at the moment of formation, the density of p-branes 
is sufficiently low (HL ^> v 2 L/£f and HL ^> cv), the p- 
brane network will experience a stretching regime. Dur- 
ing this regime, the network will be conformally stretched 
with 

L oc a and voc£f/a. (42) 

As v increases, the network will start to experience a con- 
siderable energy loss due to self-interaction (HL ~ cv) 
and, as a consequence, the Kibble regime, caracterized 
by the scaling laws 



Loz^£ f /H and v (x ^/IjH , (43) 

emerges. If, at the moment of formation of the branes, 
the density of the network is high enough, it will not 
experience the stretching regime and the Kibble regime 
may occur right away. The friction lengthscale, if, grows 



VI. CONCLUSIONS 

In this paper we developed a VOS model describing 
the dynamics of p-brane networks in homogeneous and 
isotropic backgrounds. We used it to determine the evo- 
lution of the networks' characteristic length and RMS 
velocity in flat expanding or collapsing FRW Universes, 
obtaining the corresponding scaling laws in frictionless 
and friction dominated regimes. Our model provides a 
unified semi-analitic description of p-brane network dy- 
namics, highlighting the common and distinct features 
characterizing the evolution of p-brane networks of dif- 
ferent dimensionality The connection between cosmic 
string and domain wall network evolution in 3 + 1 dimen- 
sions, which is clearly established in this work, is of par- 
ticular relevance. Although we have assumed throughout 
the paper that all p-branes have the same tension, our 
model may be extended to incorporate the possibility of 
multi-tension networks. 
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